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摘 要 : 为 研究 残 抓 沟 密 区 苹果 园 不 同 果农 复合 种 植 模 式 下 土壤 入 渗 能 力 及 其 影响 因素 ,采用 野外 双环 入 渗 法 测 
定 了 永寿 县 果 - 芒 (CK) 果 - 草 (M1) 果 - 蔬 (M2) 果 - 粮 (M3) 复 合 种 植 模式 下 全 果园 的 土壤 水 分 人 渗 参数 ,使 用 相关 性 
分 析 和 线性 匈 余 分 析 研 究 了 土壤 入 渗 特征 参数 的 影响 因素 。 利 用 3 种 常用 的 土壤 入 渗 模型 对 4 种 果园 入 渗 过 程 进 
行 拟 合 。 结 果 表 明 :4 种 果农 复合 种 植 模 式 下 果园 土壤 初始 入 渗 速 率 在 7.18~13.57 mm min '; 稳 定 入 渗 速 率 在 
1.87~2.94 mm min '; 平 均 入 渗 速 率 在 3.36~5.65 mm: min; 累计 入 渗 量 在 260.51~423.65 mm ,各 果农 复合 种 植 模 式 
苹果 园 初 始 入 渗 速 率 和 平均 入 渗 速 率 表现 为 M1>M3>M2>CK ;稳定 入 渗 速率 和 累计 入 渗 量 表现 为 M1>M2>M3>CK， 
M1 样 地 各 入 渗 指标 显著 大 于 其 他 管理 模式 样 地 ,CK 样 地 各 入 渗 特 征 指标 最 小 。 相 关 性 分 析 发 现 , 土 壤 入 渗 特 征 指 


标 与 土壤 容重 、 黏 粒 含量 和 粉 粒 含量 呈 负 相关 ;与 最 大 持 水 .最 小 持 水 量 ,毛管 持 水 量 .毛管 孔 辽 度 , 非 毛管 孔 陈 度 


总 孔隙 沙 粒 含量 呈正 相关 。 线 性 元 余 分 析 结 果 表 明 ,土壤 粉 粒 含量 、 土 壤 孔 际 度 .土壤 持 水 量 .土壤 容重 是 影响 土 
壤 人 渗 能 力 的 主要 的 因素 。 主 成 分 分 析 评价 的 土壤 和 人 渗 能 力 排序 为 M1(2.75)>M2(0.04)>M3(-0.63)>CK(-2.17) ， 


苹果 园林 下 种 植 作 物 能 够 有 效 的 增加 土壤 人 渗 能 力 。 上 | 


CASE 3 种 模型 发 现 , 欧 定 生 模 型 对 残 扰 沟 塞 区 4 种 果农 复合 


种 植 模式 果园 土壤 入 渗 过 程 的 拟 合 精度 较 高 (R=0.97) ,适用 于 描述 该 地 区 土壤 入 渗 的 实际 情况 。 
关键 词 : REAR; 复合 种 植 模式 ;土壤 入 渗 ; 人 渗 模 型 ;影响 因素 ; 陕西 


土壤 入 渗 是 指 降水 或 灌溉 水 通过 地 表 进 入 土 
壤 的 过 程 " ,是 降雨、 地 表 水 土壤 水 和 地 下 水 相互 
转化 的 重要 环节 之 一 ”。 十 壤 入 渗 直 接 影 响 着 植物 
对 土壤 水 分 的 吸收 利用 土壤 水 肥 迁 移 、 地 表 产 流 
Pow .降雨 再 分 配 等 问题 ”目前 ,国内 外 学 者 在 土 
二 入 渗 方 面 进 行 了 很 多 研究 ,并 取得 重要 成 果 , 其 
中 Kostiakov 模型 Horton 模型 .将 定 生 模型 都 是 国 
内 外 学 者 在 研究 过 程 中 提出 的 入 渗 方 程 ,对 土壤 入 
渗 过 程 均 能 够 较 好 的 阐释 。 同 时 ,关于 土壤 入 滩 方 
面 的 研究 主要 集中 于 林地 四、 草地 四 及 农田 等 自然 
生态 系统 和 人 工 生 态 系统 ,果园 作为 受 人 为 干预 的 
半 自 然 生态 系统 ,关于 土壤 人 渗 的 研究 较 少 " ,同时 
果园 内 多 种 果农 复合 种 植 模 式 对 果园 内 土壤 理化 


收 稿 日 期 : 2020-01-10; 修订 日 期 : 2020-04-27 


性 质 产 生 直接 的 影响 ”, 合 理 的 果农 复合 种 植 模式 
能 够 以 物理 手段 提升 土壤 中 质量 ,对 改善 果树 水 肥 
条 件 和 生长 环境 具有 重要 意义 ”, 而 苹果 园 果 农 复合 
种 植 模式 对 土壤 入 渗 能 力 产 生 的 影响 鲜 见报 道 。 
陕西 残 娠 沟 密 区 地 形 破 雄 , 气 候 干 燥 , 降 水 分 
布 不 均 ,是 我 国生 态 环境 最 为 脆弱 的 地 区 之 一 。 
近 几 年 , 随 着 国家 扶贫 政策 的 推进 ,苹果 产业 已 成 
为 该 地 区 主导 产业 之 一 ,在 带动 区 域 农业 经 济 中 发 
挥 着 重要 的 作用 。 然 而 ,果园 作为 耗 水 强度 和 耗 水 
深度 较 高 的 高 耗 水 作物 ,其 大 面积 种 植 加 剧 了 该 地 
区 的 水 资源 供需 矛盾 ,使 得 水 分 成 为 制约 苹果 生长 
及 产量 等 的 重要 因素 "”。 已 有 研究 表明 , HHK 
高 龄 苹果 园 土壤 训 面 普遍 存在 干燥 化 现象 , 干 层 深 
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度 其 至 可 达 10 m 以 下 。 干 层 的 存在 一 定 程度 上 
阻碍 了 雨水 对 深层 土壤 水 的 补给 ,对 苹果 生产 具 
AREER VE. FE BLT AEX AS Te] RAR EZ A 
植 模式 下 土壤 入 渗 特 征 及 其 主要 影响 因子 可 以 更 
好 地 了 解 当 前 果园 管理 模式 的 合理 性 ,为 残 抓 沟 军 
区 苹果 园 果农 复合 种 植 模式 的 选择 提供 数据 支持 ， 
对 于 苹果 园 的 可 持续 性 发 展 具 有 重要 意义 。 


1 材料 与 方法 


1.1 研究 区 概况 

实验 样 地 位 于 陕西 省 咸阳 市 永寿 县 马 坊 扰 ”， 
地 理 位 置 为 34"29'~34?59'N , 107°56’~108°20'E, 1% 
区 位 于 渭 北 高 原 南 缘 , 属 暧 温带 大 陆 性 气候 。 年 平 
均 气 温 10.2 %C ,平均 降雨 量 610.66 mm, LW 210 
d, 土 壤 母 质 为 第 四 纪 风 积 黄土 ,属于 宰 土 土 类 ,十 
土 亚 类 。 地 下 水 埋 深 20~120 m。 属 典型 的 旱 作 农 
Mb BX, ik Sor Hb JA] DCTS Se Ue. ER (Malus pumila 
Mill. ) 经 济 林 种 植 面积 较 大 , 果 业 收入 占 农民 收入 
40% 以 上 。 该 地 区 耕 层 土壤 有 机 质 含量 为 11.66 g- 
kg ', 碱 解 氮 含量 为 26.69 mg kg ,速效 磷 含 量 为 
16.94 mg 'kg ,速效 钾 含 量 为 100.59 mg*kg', pH 为 
7.65 ,容重 1.30 g.cm >。 
1.2 研究 方法 
1.2.1 样 地 选择 ”根据 果农 复合 种 植 模式 的 不 同 ， 
H T R-i (CK) (苹果 + 人 工 除草 )、 果 - 草 (M1) 
[R+ (Chenopodium album L.) + = J (Digitaria 
sanguinalis Scop.) | .5&-#i (M2) | 3258+ HK (Brassi- 
ca oleracea L.) | 、 果 - 粮 (M3)[ 苹 果 + 玉 米 (Zeu mays 
L.)j4 种 具有 代表 性 的 果农 复合 种 植 模式 果园 ,试验 


1.2.2 土壤 入 渗 能 力 测定 土壤 入 渗 试 验 于 2018 
年 7 月 在 陕西 省 咸阳 市 永寿 县 马 坊 镇 果园 进行 ,各 
AR be] 4) E |r] — UR TH 34.776° ~34.781°N, 108.162° ~ 
108.299°E) ,海拔 在 1022~1031 m, 受 环境 差异 性 影 
响 较 小 。 分 别 在 上 述 4 种 果农 复合 种 植 模式 苹果 园 
各 选取 2 个 面积 为 400 m?(20 mx20 m) 的 样 地 ,共计 
8 个 样 地 。 每 个 样 地 选取 3 个 较为 平整 的 试验 地 进 
行 双 环 和 人 渗 试 验 。 试 验 所 用 双环 内 环 直 径 10 cm, 
外 环 直 径 20 cm ,内 外 环 高 均 为 20 ecm。 试验 前 去 除 
地 表 作物 ,将 圆 环 用 橡胶 锤 缓慢 均匀 地 打 和 人 土 中 15 
cm, 尽 量 保持 土壤 结构 不 受 破坏 。 然 后 向 内 外 环 中 
同时 注水 , 当 内 环 和 外 环 水 深 达 到 3 cm 时 开始 用 秒 
表 计 时 。 并 分 别 在 0、0.5、1、2、3、5、7、10、15、20、25、 
30 min 及 之 后 每 10 min 读 取 供 水 桶 标尺 刻度 ,直到 
达到 稳 渗 状态 停止 试验 。 试 验 过 程 中 保持 内 外 环 
水 面 高 度 为 3 cm ,读数 时 一 并 读 取水 温 。 为 了 减少 
土壤 含水 量 对 试验 结果 的 影响 ,试验 时 间 相 对 集 
中 ,同时 确保 在 入 渗 实验 前 15 d 内 研究 区 未 发 生 降 
雨 。 在 土壤 入 渗 能 力 测定 同时 每 个 试验 地 按照 10 
cm 为 1 个 土屋 ,使 用 环 刀 逐 层 采集 0~50 cm + HERE 
品 各 3 个 用 于 土壤 容重 、 最 大 持 水 量 、 最 小 持 水 量 、 
毛管 持 水 量 、 非 毛管 孔 际 度 、 毛 管 孔 际 度 、 总 孔 际 度 
等 土壤 指标 的 测定 。 具 体 分 析 方 法 参照 《森林 土壤 
分 析 方 法 六 3。 各 样 地 基本 情况 见 表 1。 

1.2.3 土壤 水 分 入 渗 模 型 的 选择 ”十 壤 入 渗 模 型 是 
人 研究 土壤 入 渗 的 重要 方法 ,国内 外 学 者 建立 了 许多 
土壤 入 渗 经 验 模型 ,用 来 模拟 土壤 入 渗 随 时 间 变 化 
特征 ,本 研究 选用 3 种 比较 经 典 的 土壤 入 渗 经 验 模 
型 对 实测 数据 进行 拟 合 , 即 Horton 模型 Kostiakov 
模型 和 将 定 生 模型 ,各 模型 公式 及 参数 见 表 2。 


果园 均 位 于 同一 地 面 上 的 连 片 果园 ,品种 均 为 红 宙 
E, ERKAT (Malus prunifolia Borkh. ) ,果树 树龄 
为 20 a, 栽 植 密度 4 mx4 m, 各 复合 种 植 模式 化 肥 施 
用 量 ` 有 机 肥 施 用 量 等 农事 管理 一 致 。 


1.3 数据 统计 与 分 析 

试验 数据 使 用 Origin 2018 、SPSS 20.0 和 Canoco 
5.0 进 行 绘图 与 分 析 。 所 有 数据 均 经 正 态 分 布 检验 ， 
对 非 正 态 分 布 的 数据 进行 平方 根 变 换 。 采 用 单 因 


表 1 试验 区 土壤 理化 性 质 


Tab.1 Basic physical and chemical properties of the soil in the study area 


有 机 质 / 容重 / 


最 大 持 水 最 小 持 水 毛管 持 水 毛管 孔隙 非 毛管 孔 AILA 


模式 TEM, (g-kg"') (g-em*) 量 /mm 量 /mm 量 /mm 度 /% WE% KE% TR% WEU% DR% 
Fic CK 10.93 7.44 1.24 41.90 30.03 33.55 41.29 10.40 10.75 0.08 96.02 3.90 
果 - 草 M1 12.40 7.82 1.42 31.83 24.12 27.11 38.13 6.62 12.09 0.14 89.66 10.20 
果 - 蔬 M2 11.88 7.92 1.31 34.08 26.32 29.79 38.97 5.63 11.92 0.34 91.22 8.45 
果 - 粮 M3 12.64 7.33 1.28 38.43 29.14 32.89 41.93 TAL 11.51 0.27 90.66 9.08 
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表 2 土壤 水 分 入 渗 模 型 公式 及 参数 说 明 
Tab.2 Formulas of soil water infiltration model and 


explanation of its parameters 


模型 名 称 公式 参数 说 明 

Horton SOSSAR- 大和 大 分 别 为 初始 人 渗 速 率 
和 稳定 入 渗 速率 ;为 经 验 参 
数 ;1 为 入 渗 时 间 

Kostiakov f(D)=ar™ a 为 初始 入 渗 速 率 ;t 为 人 渗 时 
间 ;3 为 衰 碱 系数 

HE SOSLAIA /和 /分 别 为 初始 入 渗 速 率 和 


稳定 和 信 渗 速率 ;5 为 指数 ;1 为 
时 间 


素 方差 分 析 (one-way ANOVA) 对 不 同 果农 复合 种 植 
模式 果园 人 渗 参 数 进行 差异 性 分 析 。 用 了 Pearson 相 
关 性 分 析 和 线性 宛 余 分 析 (RDA) 来 表示 土壤 入 渗 特 
征 与 土壤 物理 性 质 之 间 的 关系 。 


2 结果 与 分 析 


21 不 同 果 农 复合 种 植 模式 果园 土壤 入 渗 特 征 

由 图 1 可 知 , 研 究 区 各 果农 复合 种 植 模 式 果 园 
土壤 入 渗 速 率 前 10 min 内 急速 下 降 , 人 处 于 瞬 变 阶 
段 。 此 期 间 水 分 主要 是 在 分 子 力 的 作用 下 向 下 运 
动 。 在 人 渗 开 始 10 min 时 各 样 地 平均 入 渗 速 率 较 
1 min 时 下 降 了 66%。 在 10~60 min 期 间 处 于 渗 漏 阶 


a — CK 
一 —e— M1 
g 20 —4— M2 
& 4 —— M3 
% 
网 
ES 
< 

0 T 


0 20 40 60 80 100 
时 间 /min 


累计 入 渗 量 /mm 
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段 ,入 渗 速 率 下 降幅 度 明 显 变 缓 。 同 时 ,在 此 期 间 
和信 渗 速 率 会 存在 明显 的 波动 。 人 渗 开 始 60 min 后 
土壤 入 渗 速 率 达 到 相对 稳定 的 状态 。 结 合 累 计 和 人 
渗 量 变化 来 看 ,由 于 入 渗 前 期 土壤 入 渗 速 率 较 快 ， 
因此 入 渗 开 始 时 累计 入 渗 量 增 张 幅度 较 大 ,而 到 达 
60 min 后 土壤 入 渗 速 率 达 到 稳定 状态 ,累计 入 渗 量 
斜率 不 再 变化 , 且 不 同 果农 复合 种 植 模式 土壤 累积 
入 渗 量 曲 线 斜率 差异 明显 。 

从 表 3 可 以 看 出 ,M1 样 地 各 入 渗 特 征 指标 都 是 
最 大 的 ,该 样 地 入 渗 能 力 要 高 于 其 他 样 地 。 各 果农 
复合 种 植 模式 果园 初始 入 渗 速 率 和 平均 入 渗 速 率 表 
现 为 :M1>M3>M2>CK; 稳 定 人 渗 速 率 和 累计 入 渗 量 
表现 为 :M1>M2>M3>CK, 其 中 CK 样 地 各 入 渗 特 征 
指标 与 其 他 3 种 果农 复合 种 植 模式 样 地 差异 显著 。 
2.2 土壤 入 渗 特 征 影响 因素 分 析 

土壤 入 渗 是 一 个 复杂 过 程 , 受 多 种 因素 的 影 
响 ,选取 土壤 容重 、 最 大 持 水 、 最 小 持 水 量 .毛管 持 
水 量 BFL BRE ARAB TLE .总 孔 险 度 Ata, 
含量 . 粉 粒 含量 . 沙 粒 含量 10 个 土壤 环境 因子 作为 
影响 因子 ,使 用 Pearson 相关 性 分 析 和 线性 元 余 分 析 
(RDA) 来 表示 土壤 入 渗 特 征 与 土壤 环境 因子 之 间 
KA. BRR, AHS AB PME SAE A 
ME EAA a Se ABS ; 与 其 他 土壤 性 质 呈 正 


400] — CK 

—e M1 

300] 一 M3 
250 
200 
150 
100 
50 
0 

0 20 40 60 80 100 


时 间 /min 


图 1 不 同 果农 复合 种 植 模式 土壤 入 渗 过 程 


Fig. 1 Soil infiltration process under different fruit—agriculture intercropping patterns 


表 3 不 同 果农 复合 种 植 模 式 土壤 入 渗 特 征 


Tab.3 Soil infiltration characteristics under different fruit-agriculture intercropping patterns 


复合 种 植 模式 初始 入 渗 速率 /mm.min) ”稳定 入 渗 速率 /(mm:min"') ”平均 入 渗 速率 /(mmmin') 累计 入 渗 量 /mm 
CK 7.18+0.09c 1.87+0.02c 3.36+0.07c 260.51+4.26d 
MI1 13.5740.29a 2.94+0.09a 5.65+0.12a 423.65+6.32a 
M2 11.29+0.16b 2.34+0.06b 4.56+0.14b 335.71+8.82b 
M3 13.38+40.14a 2.15+0.09b 4.74+0.13b 320.58+5.28c 


注 : 同 列 数 据 后 不 同 小 写字 母 表示 差异 显著 (P<0.05 )。 
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表 4 土壤 入 渗 特 征 与 土壤 因子 相关 关系 和 矩阵 


Tab.4 Matrix of correlation between soil infiltration characteristics and soil factors 


指标 容重 ”最 大 持 水 量 最 小 持 水 量 毛管 持 水 量 毛管 孔 际 度 ” 非 毛 管 孔 际 度 AIRE Abe 粉 粒 含量 沙 粒 含量 
初 渗 速 率 -0.886" 0.821” 0.843” 0.901" 0.847" 0.385 0.774  -0.117 -0.890” 0.878” 
稳 渗 速率 -0.596 0.664° 0.562 0.495 0.337 0.707 0.552 -0.345 -0.457 0.443 
平均 入 渗 速 率 -0.799” 0.805™ 0.762° 0.757 0.645" 0.543 0.724 -0.133  -0.756" 0.743” 
累计 和 人 渗 量 -0.668- 0.720" 0.634" 0.585 0.44 0.6 0.617 -0.289 -0.805" 0.801” 


注 :* 表 示 在 0.05 水平 上 显著 相关 ;** 表 示 在 0.01 水平 显著 相关 。 


BELE 初始 人 


图 2 土壤 入 渗 特征 与 土壤 因子 RDA 排序 
Fig. 2 Soil infiltration characteristics and 


soil factor RDA ranking map 


相关 ( 表 4, 图 2)。 
从 Pearson 相关 性 分 析 结 果 来 看 ( 表 4) ,土壤 初 
渗 速 率 .平均 人 活 速 率 与 土壤 容重 和 粉 粒 含量 呈 极 


表 5 土壤 因子 对 土壤 入 渗 特 征 影响 贡献 度 


Tab.5 Contribution of soil factors to soil infiltration 


characteristics 

土壤 因子 解释 率 /% F P 
粉 粒 含量 54.8 18.4 0.002™ 
总 孔隙 度 18.2 3.4 0.084 
非 毛 管 孔 际 度 7.7 2.5 0.126 
最 小 持 水 量 6.7 5.7 0.044" 
最 大 持 水 量 5.6 2.9 0.116 
AE E 4.3 7.8 0.046" 
容重 1.4 26.1 0.02* 
毛管 孔 际 度 1.2 3 0.166 
沙 粒 含量 0.1 8.4 0.194 


BE .最 小 持 水 量 .最 大 持 水 量 。 说 明 土 壤 颗 粒 组 E 
RILE .土壤 持 水 能 力 是 影响 土壤 人 渗 特 征 的 主 
23 不 同 果农 复合 种 植 模式 果园 土壤 入 渗 能 力 
评价 

为 了 评价 不 同 果农 复合 种 植 模式 下 残 坟 沟 窜 
区 苹果 网 土壤 入 渗 能 力 ,采用 SPSS 软件 对 初始 入 渗 


显著 负 相 关 ; 累计 人 渗 量 与 容重 呈 显 著 负 相关 ,与 
粉 粒 含 量 呈 极 显 著 负 相关 。 容 重 越 大 , 粉 粒 含量 越 
多 ,入 渗 性 能 越 差 。 入 渗 特 征 指 标 与 最 大 持 水 量 、 
最 小 持 水 量 .毛管 持 水 量 .毛管 孔隙 度 , 非 毛 管 孔 院 
度 .总 孔隙 度 ` 沙 粒 含量 呈正 相关 。 黏 粒 含量 对 该 
地 区 土壤 入 渗 特 征 影响 均 不 显著 。 

为 进一步 说 明 究竟 哪 种 土壤 因子 与 人 渗 特 征 
指标 关系 更 密切 ,采用 线性 元 余 分 析 对 10 个 土壤 因 
子 与 人 渗 关 系 进行 分 析 。 从 土壤 因子 影响 土壤 入 
渗 特 征 的 RDA 分 析 结 果 来 看 ( 表 5, 图 2) ,本 文 所 选 
指标 共 人 解释 土壤 入 渗 特 征 的 99.4%。 其 中 土壤 粉 粒 
含量 与 各 入 渗 指 标 负 相关 性 最 强 ;土壤 非 毛 管 孔 际 
度 与 各 入 渗 指 标 正 相关 性 最 强 ( 图 2)。 各 土壤 环境 
因子 对 土壤 入 渗 特 征 影响 贡献 度 ( 表 5) 表 现 为 粉 粒 
含量 最 大 ,为 54.8%, 其 次 为 总 孔 际 度 、 非 毛管 孔隙 


速率 (XX) FARE A BIR (XD) .平均 入 渗 速 率 (X)、 
累计 入 渗 量 ( 生 )4 个 因素 进行 主 成 分 分 析 , 得 出 主 
成 分 方程 ,计算 各 复合 种 植 土壤 渗透 能 力 得 分 。 表 
6 显示 2 个 主 成 分 贡献 率 达 到 98.28% ,可 以 准确 的 
对 总 方差 进行 解释 ,能 够 解释 为 土壤 入 渗 能 力 的 量 
度 。 其 主 成 分 方程 为 P,=0.491X1+0.458+0.521Xs+ 
0.528X4; Ps=—0.545X,+0.740X,-0.340X;+0.200X。 
将 标准 化 后 的 土壤 入 渗 特 征 参 数 代 入 主 成 分 
方程 ,计算 得 到 不 同 果农 复合 种 植 模 式 下 各 样 地 的 


表 6 土壤 水 分 入 渗 能 力 的 PCA 分 析 
Tab.6 PCA analysis of soil water infiltration capacity 
累计 贡献 
率 /% 
P, 0.491 0.458 0.521 0.528 3.444 86.09 86.09 
P, -0.545 0.740 -0.340 0.200 0.487 12.19 98.28 


主 成 分 X X% Xs X 特征 值 贡献 率 /% 
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表 7 土壤 入 渗 能 力 评价 


Tab.7 Evaluation of soil infiltration capacity 


复合 种 植 模式 P, P, 得 分 排序 
CK -2.41 0.24 -2.17 4 
M1 2.28 0.47 2.75 1 
M2 0 0.04 0.04 2 
M3 0.13 -0.76 -0.63 3 


得 分 , 表 7 显 示 , 不 同 果 农 复 合 种 植 模式 下 土壤 入 渗 
能 力 强 弱 排序 为 : 果 - 草 (M1)> 果 - 蔬 (M2)> 果 - 粮 
(M3)> 果 -项 (CCK)。 果 - 草 (MI)、 果 - 蔬 (M2)、 果 - 粮 
(M3) 样 地 因 存 在 草本 ,使 得 植物 根系 较 多 ,地 表土 
质 芯 松 ,所 以 土壤 人 渗 能 力 较 强 。 果 - 荒 (CCK) 样 地 
由 于 除草 的 原因 ,植物 根系 稀少 ,土壤 结构 差 , 所 以 
土壤 入 渗 能 力 最 弱 。 
2.4 入 渗 速 率 的 模型 拟 合 与 评价 

由 于 土壤 人 渗 是 在 一 定时 间 内 逐渐 减缓 的 过 
程 ,土壤 的 初始 入 渗 速 率 很 大 ,而 稳定 入 渗 速 率 则 
较 小 。 不 同 果 农 复 合 种 植 模式 下 ,不 同 入 渗 模 型 对 
土壤 的 初始 入 渗 和 稳定 入 渗 阶 段 拟 合 的 效果 也 不 
相同 ,因此 ,本 研究 使 用 3 种 常用 的 土壤 入 渗 经 验 模 
型 对 实测 数据 进行 拟 合 ( 表 8)。 


重 凡 旭 等 :不 同 果农 复合 种 植 模式 土壤 人 渗 能 力 及 其 影响 因素 
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表 8 显示 ,不 同 果农 复合 种 植 模式 土壤 入 渗 过 
程 的 模型 拟 合 优 度 存 在 差异 ,从 大 到 小 为 :将 定 生 
模型 >Kostiakov 模型 >Horton 模型。 其 中 ,将 定 生 模 
型 的 决定 系数 为 0.96 ~ 0.99, 平 均值 为 0.97, Kostia- 
kov 模型 的 决定 系数 为 0.93 ~ 096 ,平均 值 为 0.95。 
Horton 模型 的 决定 系数 为 0.91 ~ 0.96 ,平均 值 为 
0.93。 将 定 生 模型 能 较 好 的 反应 不 同 果农 复合 种 植 
模式 下 的 土壤 水 入 渗 的 实际 情况 ,可 以 更 准确 的 对 
残 卉 沟 密 区 不 同 果 农 复合 种 植 模式 下 苹果 园 土 壤 
水 的 入 渗 过 程 和 入 渗 能 力 进 行 模拟 分 析 。 

使 用 平均 相对 误差 和 均 方 根 误差 对 各 模型 拟 
合 优 度 进行 对 比分 析 ( 表 9) ,结果 显示 :Horton 模 
型 Kostiakov 模型 . 蒋 定 生 模 型 的 平均 相对 误差 的 
变化 范围 分 别 是 5.49% ~ 39.95% .7.93% ~ 12.36% 、 
7.31% ~ 15.33% , 均 方 根 误差 的 变化 范围 分 别 是 
0.38~1.42 mm.min ',0.40~1.01 mm: min ,0.26~0.87 
mm min!。 对 比 3 个 模型 可 知 , 蒋 定 生 模型 的 相对 
误差 较 小 ,与 土壤 入 渗 的 实测 过 程 更 贴近 。 因 此 ， 
结合 各 个 样 地 土壤 入 渗 模 型 拟 合 优 度 和 平均 相对 
误差 的 大 小 ,认为 蒋 定 生 模型 更 符合 残 壕沟 殉 区 土 
培 入 渗 的 实际 过 程 ,可 以 更 好 地 模拟 其 和 人 渗 过 程 。 


HS 土壤 入 渗 模 型 拟 合 参 数 


Tab.8 Fitting parameters of soil infiltration model 


Horton 模型 Kostiakov 模 型 蒋 定 生 模 型 
复合 种 植 模式 
fo k R a b R J: fi b R 
CK 1.90 7.84 0.08 0.96 7.69 0.30 0.96 =5.97 7.65 0.10 0.98 
M1 4.02 26.95 0.62 0.91 15.29 0.39 0.95 219 12.87 0.57 0.96 
M2 3.29 26 0.82 0.91 12.93 0.42 0.93 2.27 10.18 0.72 0.96 
M3 2.84 22.50 0.41 0.92 14.68 0.46 0.95 1.04 12.44 0.55 0.99 
表 9 实测 入 渗 过 程 与 模拟 值 误差 比较 
Tab.9 Error comparison of the measured soil infiltration rate and the calculated value 
Horton 模型 Kostiakov 模型 FEAL 
复合 种 植 模式 
MRE/% RMSE/(mm- min" ) MRE/% RMSIE/(mm: min) MRE/% RMSIE/(mm- min") 
CK 5.49 0.38 12.24 0.40 9.07 0.26 
M1 17.19 1.39 10.51 1.01 9.40 0.87 
M2 39.95 1.42 12.36 1.01 15.33 0.79 
M3 25.28 1.22 7.93 0.51 Fol 0.43 


HE: MRE 为 平均 相对 误差 ;RMSE 为 均 方 根 误差 。 
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ar 草 的 密度 。 同 时 农作物 和 蔬菜 会 在 其 成 熟 后 移出 
3 讨论 果园 , 仅 剩余 植物 的 根系 。 因 此 果蔬 间作 与 果农 间 


土壤 水 分 人 渗 过 程 是 水 分 在 分 子 力 .毛管 力 和 
重力 的 综合 作用 下 在 土壤 中 运动 的 物理 过 程 ,入 渗 
按 土壤 水 分 受 力 情况 及 运行 特征 可 分 3 个 阶段 "， 
即 渗 润 阶段 渗 漏 阶 段 和 渗透 阶段 段 。 本 文 研究 发 
现 , 不 同 果农 复合 种 植 模式 果园 土壤 入 渗 过 程 存在 
差异 ,但 是 都 具有 初 渗 速 率 > 平均 人 渗 速 率 > 稳定 
和信 渗 速率 。 究 其 原因 , 渗 润 阶段 水 分 主要 受 分 子 力 
作用 ,被 土壤 颗粒 吸附 而 形成 薄膜 水 ,对 于 干燥 的 
土壤 ,该 阶段 入 渗 率 明显 大 于 其 他 阶段 。 渗 漏 阶段 
下 渗水 主要 受到 重力 和 毛细 管 引力 ,该 阶段 土壤 水 
主要 在 土壤 颗粒 间 的 孔隙 移动 ,并 开始 充填 孔 际 ， 
直至 充满 水 ,此 阶段 下 渗 速 率 迅速 递减 。 渗 透 阶 段 
土壤 孔隙 已 被 水 充满 达到 饱和 ,水 分 在 重力 作用 下 
移动 ,属于 饱和 水 流 运动 ,水 分 人 渗 率 维持 稳定 。 

对 果园 入 渗 方面 的 研究 大 多 是 从 不 同类 型 果 
园 方 面 进行 ” ,但 是 对 于 不 同 果农 复合 种 植 模 式 下 
果园 入 渗 能 力 的 变化 报道 较 少 。 农 林 生 产 过 程 中 
多 样 的 管理 方式 具有 改善 生态 系统 服务 功能 的 作 
用 外 ,同时 也 会 改变 土壤 理化 性 质 ,直接 或 间接 影响 
着 土壤 入 渗 能 力 。 苹 果园 林 下 作物 及 其 根系 分 布 
对 土壤 粒度 组 成 . 筷 际 度 等 物理 性 质 有 着 较 大 的 改 
善 作用 2 , 林 下 作物 不 仅 可 以 在 一 定 程度 上 削减 
雨滴 对 土壤 的 打击 力 ,还 能 够 分 解 为 土壤 所 需要 的 
养分 ,改善 土壤 结构 ;较为 发 达 的 植物 根系 可 以 
将 有 机 质 纳 入 土壤 ,同样 也 能 够 改善 土壤 结构 。 
根系 在 死亡 腐烂 后 会 使 土壤 存在 优先 通道 3, 进 而 
增加 土壤 水 流 的 下 渗 , 对 入 渗 起 到 促进 作用 。 果 农 
复合 种 植 模式 的 变化 同样 影响 着 小 型 土壤 动物 的 
变化 , 林 下 生 草 可 以 明显 的 增加 小 型 土壤 动物 的 数 
量 汪 ,小 型 土壤 动物 的 存在 对 土壤 性 质 的 改善 同样 
也 起 到 了 一 定 的 作用 。 而 CK 样 地 在 果树 生长 期 间 
在 除草 时 会 对 土壤 表面 践踏 ,这 也 是 土壤 入 渗 速 率 
降低 的 主要 原因 。 本 文 研究 还 发 现 , 林 下 存在 草本 
层 的 各 样 地 土壤 入 渗 能 力 明显 要 优 于 除草 样 地 ,其 
结果 可 归结 于 不 同 果农 复合 种 植 模式 下 草本 根系 
生物 量 、 根 系 活动 .土壤 动物 的 变化 以 及 人 为 活动 
影响 而 导致 。 同 时 M2 和 M3 样 地 土壤 入 渗 速 率 也 
小 于 M1l 样 地 ,其 原因 主要 是 虽然 M2 和 M3 林 下 存 
在 植被 ,但 是 M2 和 M3 林 下 植被 的 密度 要 远 小 于 生 


作 样 地 土壤 入 渗 能 力 要 小 于 生 草 样 地 。 

不 同 果农 复合 种 植 模式 果园 ,土壤 理化 特性 也 
存在 相应 差异 ,直接 或 间接 影响 着 土壤 入 渗 能 力 。 
本 文 使 用 相关 性 分 析 和 RDA 分 析 对 影响 入 渗 特 征 
的 土壤 物理 性 质 进行 了 分 析 ,结果 表明 土壤 颗粒 组 
成 。 了 和 孔隙 度 和 容重 是 影响 该 地 区 土壤 入 渗 特 征 的 
儿 个 最 主要 因素 。 其 中 粉 粒 含量 容重 AMA E 
对 土壤 入 渗 起 着 负 向 的 影响 , 即 粉 粒 含 量 BL 
量 和 容重 越 大 土壤 人 渗 能 力 越 差 ;其 他 几 个 指标 对 
土壤 入 渗 起 着 正身 的 影响 。 这 与 吕 刚 等 守 和 李 蛙 29 
的 结果 一 致 。 在 各 指标 中 和 粒 含量 对 土壤 入 活 的 
负 癌 影响 不 显著 ,其 原因 主要 是 由 于 人 研究 区 土壤 以 
粉 粒 含 量 为 主 ,其 次 为 沙 粒 含 量 , 笑 粒 含 量 所 占 比 
例 较 小 ,含量 仅 在 0.08%~0.34% 之 间 ( 表 1)。 

土壤 入 渗 模 型 是 研究 土壤 入 渗 的 重要 方法 , 土 
二 入 渗 模 型 能 够 表达 入 渗 速 率 随时 间 变 化 规律 ,从 
而 对 土壤 入 渗 过 程 作出 定性 或 定量 的 评价 。3 种 常 
HERA IT Gi Vad AX DK SE ER ro] BE AB BI BE 
REMAIN A FE rE AE BAY. Kostiakov 模型 Horton 
模型 ,其 中 Kostiakov 模型 是 假设 起 始 时 的 人 渗 速 率 
是 无 穷 大 2 ,当时 间 无 限 增 大 时 ,入 渗 速 率 将 趋 近 
于 0。 但 在 本 文中 ,由 于 在 垂直 和 人 渗 的 情况 下 有 重 
力 势 的 作用 ,入 渗 速 率 最 后 会 接近 于 稳定 人 渗 率 后 
保持 稳定 。 而 蒋 定 生 模型 和 Horton 模型 均 增 加 了 
常数 项 以 反映 无 限 长 时 在 重力 作用 下 达到 了 稳定 
和信 渗 下 ,这 更 符合 土壤 水 动力 学 。 将 定 生 模型 和 
Horton 模型 2 个 模型 中 , 蒋 定 生 模 型 具有 更 大 的 拟 
合 优 度 , 因 此 ,将 定 生 模型 更 适宜 描述 残 姑 沟 密 区 
苹果 园 土壤 入 渗 过 程 ,这 与 吴 秋 菊 等 ”在 黄土 区 的 
研究 结果 一 致 。 开 展 不 同 果 农 复 合 种 植 模 式 对 土 


培 入 渗 能 力 变化 的 研究 可 以 为 果园 管理 方式 的 选 
择 提供 依据 。 
4 结论 


通过 野外 双环 入 渗 实 验 对 残 杖 沟 罕 区 苹果 园 4 
种 果农 复合 种 植 模式 土壤 入 渗 过 程 进行 观测 ,评价 
3 种 常用 入 渗 模 型 在 残 堪 沟 蜜 区 果园 适用 性 ,分 析 
容重 等 10 个 土壤 环境 因子 对 土壤 入 渗 特征 指标 的 
影响 ,最 后 使 用 主 成 分 分 析 法 评价 了 4 种 苹果 园 土 
二 入 渗 能 力 。 通 过 分 析 得 出 以 下 4 个 结论 : 
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(1) 4 种 果农 复合 种 植 模式 果园 土壤 初始 人 渗 
速率 在 7.18~13.57 mm: min” 之 间 ; 稳 定 入 渗 速 率 在 
1.87~2.94 mm ' min 之 间 ; 平 均 人 渗 速 率 在 3.36~ 
5.65 mm: min’! Z [i] ; ib AB tet FE 260.51~423.65 
mm*min "之 间 , 各 果农 复合 种 植 模式 果园 和 人 渗 特 征 
参数 存在 显著 差异 。 

(2) 相关 性 分 析 表 明 :在 所 选取 的 10 个 土壤 环 
境 因 子 中 。 土 壤 人 渗 特 征 指标 与 土壤 容重 AS 
量 和 粉 粒 含量 呈 负 相关 ; 与 最 大 持 水 .最 小 持 水 量 、 
毛管 持 水 量 、 毛 管 孔 际 度 、 非 毛管 孔 际 度 、 总 孔 际 沙 
粒 含量 呈正 相关 。RDA 分 析 表 明 ;10 个 土壤 环境 因 
子 共 解释 土壤 人 渗 特 征 总 变量 的 99.4%。 各 十 壤 因 
子 中 ,土壤 粉 粒 含 量 .土壤 孔 际 度 .土壤 持 水 量 、 土 
二 容重 是 影响 该 地 区 苹果 园 土壤 入 渗 能 力 的 4 个 主 
要 的 因素 。 

(3) 主 成 分 分 析 表 明 所 选用 的 4 个 人 渗 特 征 指 
标 对 土壤 入 渗 能 力 的 贡献 量 达到 了 98.28% ,不 同 果 
农 复 合 种 植 模 式 果 园 土 壤 人 渗 能 力 排 序 :MIl 
(2.75)>M2(0.04)>M3(-0.63)>CK(-2.17) ,苹果 园 
林 下 种 植 作物 能 够 有 效 的 增加 土壤 入 渗 能 力 。 

(4) 不 同 果 农 复合 种 植 模式 土壤 入 渗 过 程 的 模 
型 拟 合 优 度 存 在 一 定 差异 ,其 拟 合 优 度 依 次 为 蒋 定 
生 模 型 (0.97) >Kostiakov 模型 (0.95) >Horton 模 型 
(0.93 ) ,其 中 蒋 定 生 模型 较 符 合 黄 土 残 震 沟 餐 区 苹 
果园 土壤 人 渗 的 实际 情况 。 
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Soil infiltration ability and influencing factors under different intercropping 
patterns in orchards in the gully region of Loess Plateau 


MENG Fan-xu'’’, WANG Shu-sen'**, MA Ying-mei””, QIN Fu-cang™, 
LUO Yu-yang'*’, WANG Di-hai’, GUO Yu', GUO Xiang-dong' 

(1. College of Desert Control Science and Engineering, Inner Mongolia Agricultural University, Hohhot 010018 
Inner Mongolia, China; 2. Key Laboratory of State Forest Administration for Desert Ecosystem Protection and 
Restoration, Hohhot 010018, Inner Mongolia, China; 3. Inner Mongolia Key Laboratory of Aeolian Physics and 
Sand Control Engineering, Hohhot 010018, Inner Mongolia, China; 4. College of Forestry, Northwest A&F 
University, Yangling 712100, Shaanxi, China) 


Abstract: We investigated the soil infiltration characteristics and influencing factors under different fruit- 
agriculture intercropping patterns below apple trees (Malus) in the gully region of Loess Plateau. A double- ring 
infiltration method was used to measure the soil water infiltration parameters under fruit- agriculture intercropping 
patterns of Malus—wild (CK), Malus—Herbage (M1), Malus—Brassica (M2), Malus—Zea mays (M3) in the orchards of 
Yongshou County. Correlation and redundancy analyses (RDA) were used to determine the influencing factors of 
soil infiltration parameters. Three common soil infiltration modells were used to fit the infiltration processes of four 
orchards. Under four fruit-agriculture intercropping patterns, the initial orchard soil infiltration rate was 7.18-13.57 
mm: min", the stable infiltration rate was 1.87-2.94 mm: min", the average infiltration rate was 3.36-5.65 mm: min", 
and the cumulative infiltration amount was 260.51—423.65 mm. The initial infiltration rate and average infiltration 
rate of the four orchards were ranked from highest to lowest as M1>M3>M2>CK; the steady infiltration rate and 
cumulative infiltration were ranked as M1>M2>M3>CK. The infiltration index of M1 was significantly higher than 
that of the other fruit- agriculture intercropping patterns, and the infiltration characteristic index of CK was the 
lowest. The correlation analysis showed that the soil infiltration characteristics were negatively correlated with soil 
bulk density, clay content, and silt content, and positively correlated with maximum water holding capacity, 
minimum water holding capacity, capillary water holding capacity, capillary porosity, non- capillary porosity, and 
total pore sand content. The RDA analysis revealed that the soil silt content, soil porosity, soil water holding 
capacity, and soil bulk density were the main factors influencing soil infiltration performance. Soil infiltration ability 
was ranked as M1 (2.75)>M2(0.04) >M3(—0.63) >CK(- 2.17) by principal component analysis. Therefore, planting 
crops under apple orchards could effectively increase soil infiltration ability. Comparing the three models, Jiang 
Dingsheng’ s model had a high fitting accuracy (R’=0.97) for the soil infiltration processes of orchards under four 
fruit- agriculture intercropping patterns in the gully region of Loess Plateau, which was suitable for describing the 
actual situation of soil infiltration in this area. 

Keywords: gully region of Loess Plateau; fruit- agriculture intercropping patterns; soil infiltration; infiltration 


model; influence factor; Shaanxi 


